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Abstract 

Current searches for doubly charged Higgs bosons (-ff ±=t ) at the Fermilab Tevatron are 
sensitive to single production of -ff ±=t , although the pair production mechanism qq — > 
H ++ H is assumed to be dominant. In the context of a Higgs Triplet Model we study 
the mechanism q'q — > H^H^ at the Tevatron and CERN Large Hadron Collider, and 
show that its inclusion can significantly improve the search potential for i/ ±=t . Moreover, 
assuming that the neutrino mass is generated solely by the triplet field Yukawa coupling 
to leptons, we compare the branching ratios of — > Z Z and — ► H^W* for 
the cases of a normal hierarchical, inverted hierarchical and degenerate neutrino mass 
spectrum. 
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1 Introduction 



The quest for Higgs bosons is of utmost importance at high energy colliders [I],(2j, jHj- In the 
Standard Model (SM), one isospin I = 1/2, hypercharge Y = 1 complex scalar doublet breaks 
the Electroweak Symmetry (EW) and provides mass for the fermions, and Z. One neutral 
scalar, 0°, remains as a physical degree of freedom - "the SM Higgs boson". Such a framework 
predicts p (= M^/M| cos 2 Ow) = 1 at tree-level, a result which is in impressive agreement with 
the experimental measurement of p ~ 1 jlj. More generally, any Higgs sector composed solely 
ofi=l/2,V = l doublets assures p = 1 at tree- level, with calculable 1-loop corrections 

Predicting p — 1 at tree-level is certainly an attractive feature of i = 1/2, Y = 1 doublet 
representations, although models with isospin triplets (i = 1) can also be considered pQ. Such 
models have various virtues and deficiencies. If the neutral member of the triplet acquires a 
vacuum expectation value (VEV) then p = 1 at tree-level is no longer guaranteed, and the 
triplet VEV must be very small in order to comply with the measured value p ~ 1. However, 
unlike doublets, Y = 2 triplets can give rise to neutrino masses and mixings whose magnitude 
is proportional to the triplet vacuum expectation value multiplied by an arbitrary Yukawa 
coupling (hij) without invoking a right handed neutrino [H], [7j. 

A clear phenomenological signature of Y = 2 triplets would be the observation of a doubly 
charged Higgs boson i/ ±± . Such if ±± have been searched for at the e + e~ collider LEP, resulting 
in mass limits of the order m H ±± > 100 GeV 0,0, [TU].|llj. Their existence can also affect 
a wide variety of processes, such as Bhabha scattering, the anomalous magnetic moment of 
the muon ((7 — 2)^, and lepton flavour violating p^ 1 and decays |12| .|13]. |14| . |15| .|16]. |17| . 
The Fermilab Tevatron recently performed the first search for if ±=t at hadron colliders. The 
production process pp — > 7, Z — > H ++ H was assumed, with subsequent decay if ±=t — > 
l ± l ± . DO [TBI searched for ii ±± -> p ± p ± while CDF PU searched for 3 final states if ±± -> 
p ± p ± , p ± e ± and e ± e ± . Mass limits of the order m H ±± > 130 GeV were obtained with an 
integrated luminosity of 240 pb , assuming BR(if ±± -> iflf) = 100% [HI in a given channel. 
These are the strongest direct mass limits on any type of Higgs boson, which shows the strong 
search capability of hadron colliders in the channel if ±=t — > l^l^. 

Given this strong search potential, in this paper we consider the phenomenological effect 
of relaxing these simplifying assumptions for the dominant production mechanism and decay 
modes of if ±± . Although work along these lines has appeared previously [H], |2"U].[2*T].[22j. 
j22], [23] we develop and expand the preceding analyses. For example, if hij are solely responsible 
for the currently favoured form of the neutrino mass matrix then BR(ii ±± -> iflf) < 100% 
in a given channel ^3] . In this paper we study in detail the alternative production mechanism 
q'q — > W* — > H^H^ (23], which can be as large as qq — > 7, Z — > H ++ H . Since the current 
search strategy at the Tevatron is in fact sensitive to single production of if ±=t , we introduce 
the inclusive single production cross-section (a^±±) as the sum of the single and pair production 
cross-sections. We point out that the contribution of q'q — > W* — * H ±:k H T to cr H ±± strengthens 
the Tevatron mass limit on ii ±± , which in general has a dependence on m#±. Moreover, we 
quantify the impact of the potentially important decay mode if ±± — > H^W* [22J in the light 
of recent neutrino data. Although such a decay can weaken the if ±=t search capability in the 
leptonic channel, observation of H ±± — > H ± W* together with one or more leptonic channels 
might permit an order of magnitude estimate of [2*3] . 



Our work is organized as follows. In Section 2 we introduce the Higgs Triplet Model. In 
Section 3 we study the production mechanism q'q — > if ±± iJ =F and its phenomenological effect 
on the H ±zL search at the Tevatron and LHC In Section 4 we quantify the impact of the 
decay H ±zL — > H ± W*, while Section 5 considers the search potential of the Tevatron in the 
generalized scenario. Finally, in Section 6 we present our conclusions. 



2 The Higgs Triplet Model 

Higgs 1 = 1 triplet representations arise in several well motivated models of physics beyond 
the SM P,pOj- For example, Left-Right (L-R) symmetric models built on the gauge group 
SU(2)r x SU(2)l x U(l) contain both left- and right-handed I = 1,Y = 2 triplet representa- 
tions. Such models also require extra gauge bosons and can provide naturally light neutrino 
masses via the seesaw mechanism. Little Higgs models [23] also require I = 1, Y = 2 triplet 
representations, as well as new gauge bosons and fermions. However, Higgs triplets can be 
considered as a minimal addition to the SM j2E|- for a review see |2Zj • We will focus on a par- 
ticularly simple model 0, [I] which merely adds a I = 1,Y = 2 complex (left-handed) Higgs 
triplet to the SM Lagrangian, hereafter referred to as the "Higgs Triplet Model" or "HTM". 
Such a model can provide a Majorana mass for the observed neutrinos without the need for a 
right handed neutrino via the gauge invariant Yukawa interaction: 1 

C = hijil/[ L CiT 2 Ail)j L + h.c (1) 

Here hij(i,j = 1,2,3) is an arbitrary coupling, C is the Dirac charge conjugation operator, 
4>iL = (viji) 1 ^ is a left-handed lepton doublet, and A is a 2 x 2 representation of the Y = 2 
complex (left-handed) triplet fields: 

A-( A+ /^ A++ (2) 
A "l A -A+/V2 {Z) 



The Higgs potential [Ej is as follows, with $ 



kO\T. 



V = m 2 ($ t $) + Ai($ t $) 2 + M 2 Tr(A t A) + A 2 [Tr(A t A)] 2 + A 3 Det(A t A) 

+A 4 ($ t $)Tr(A t A) + A 5 ($ t r l $)Tr(A t r i A) + (^=^ T ir 2 A^) + h.cj (3) 

The term /i$A$, where fi is a dimensionful trilinear coupling, gives rise to a VEV va for the 
neutral member of the triplet A : 

v A ~ /j,v 2 /2M 2 (4) 

Here M is the common triplet mass (M 2 A^A). Since we are interested in the case of light 
triplets we take M « v, and so v a ~ fi- A non-zero va gives rise to the following mass matrix 
for neutrinos: 

rriij = 2hij(A°) = V^hijVA (5) 
1 Note that the analogous term for a Y=0 triplet is forbidden by gauge invariance 
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Note that the HTM is free from a massless Goldstone boson (Majoron) arising from the violation 
of the lepton number (L) global symmetry, because the Higgs potential contains the term /x$A$ 
term which explicitly violates lepton number when A is assigned L = —2. Cosmological data 
provides a constraint on the neutrino masses TOj, Em, £0.75 eV [2H1- Lepton flavour violating 
(LFV) processes involving p, and r provide the strongest upper limits on hij and hence v a 
cannot be arbitrarily small if the HTM is to accommodate the currently favoured form of the 
neutrino mass matrix. A rough lower bound v& ;> 10 eV can be derived. An upper limit on v& 
can be obtained from considering its effect on p. In the HTM p is given by (where x = Va/v): 

From the measurement of p ~ 1 a purely tree-level analysis gives the bound v^/v £0.03. We 
will comment on the 1-loop expression for 5p below 29J, [20], [31]. in this paper we will assume 

10 eV £v A £10000 eV (7) 

Hence the tree- level value of p is essentially equal to 1, thus easily satisfying the experimental 
constraint on 5 p. Such small values of v& can be explained by a 2 loop mechanism ^3] or in 
the context of extra dimensions [22] , jUII] • Moreover, such values of v& would permit some 
to be sufficiently large to enhance various LFV p and r decays to the sensitivity of current and 
forthcoming experiments P3J 5 [To].|17j. and are also consistent with the requirement that any 
primordially generated baryon asymmetry is not erased by the lepton number violating triplet 
interactions |3"fj . 

The HTM has 7 Higgs bosons {H ++ , H~, H+, H~, H°, A , h°). While is purely triplet 
(= A ±± ), the remaining eigenstates would in general be mixtures of the doublet and triplet 
fields. Such mixing is proportional to the triplet VEV, and hence small even if v a assumes 
its largest value of a few GeV. Therefore the first six eigenstates are essentially composed of 
triplet fields, while the 7 = 1/2 doublet gives rise to a SM like h° and the Goldstone bosons 
G^jG . The most striking signature of the HTM would be the observation of i7 ±± . 2 In the 
HTM there exists the following relationships among the masses of the physical Higgs bosons: 



m\ ±± ~ M 2 + 2 [M g2 A5) M 2 w (8) 

m 2 H ± ~ m 2 H ±± + 2^-M^ 

m 2 H0 A0 ~ m 2 H ± + 2^Mw 
9 

Here M is the triplet mass term, while A4, A5 are dimensionless quartic couplings. For A5 > 
(A5 < 0) one has the following hierarchy mu±± < tuh± < fnH°,A° ( m #±± > m H ± > m H°,A )- 
Clearly M sets the scale for the mass of the triplet fields, while the mass splitting among the 

2 The dominantly triplet eigenstates H ± , H° and A can have a different phenomenology to the analogous 
Higgs bosons in doublet (/ = 1/2, Y — 1) representations. 
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eigenstates is determined by the quartic couplings and can be O(Mw)- We will focus on Higgs 
boson masses of interest for the Tevatron and LHC, and hence we assume M J> 1 TeV. 

At the 1-loop level the Higgs sector contribution to 5p is a function of v& and the Higgs 
boson masses. Although a quantitative analysis in the context of the HTM is still lacking, 
explicit formulae for the contributions of Y = 2 triplets to the self-energies of the W and Z 
in the context of L-R symmetric models and Little Higgs Models can be found in |30j. 
In particular, such contributions are sensitive to the mass splittings of the Higgs bosons. In 
the HTM the triplet Higgs boson mass splitting is determined by the quartic coupling A 5 , with 
A5 = giving rise to degenerate triplet scalars of mass M. We will present results for both the 
degenerate case and for mild splittings of up to 20 GeV in our discussion of H ±zL phenomenology 
at the Tevatron. 

We now briefly discuss present mass bounds on the Higgs bosons of the HTM, which differ 
in some cases from the commonly quoted mass bounds in the 2HDM. If H° and A were 
the lightest, they could have been produced at LEP via the mechanism e + e~ — > A°H° (note 
that e + e~ — > ZH° is proportional to v& and hence negligible). However, since A and H° 
would both decay invisibly to uV, Ref. [BE} suggested using LEP data on juV (where 7 arises 
from bremsstrahlung from e + or e~) and derived the mass limit m#o ^0 ^55 GeV. Concerning 
H ^ , LEP searched for H ± — > cs or rv T which are expected to be the dominant decays in 
doublet models, and obtained mass limits around m H ± J>80 GeV. For the triplet the 
decays — > e^v^^v may have large BRs. However, in this scenario one could presumably 
use data from slepton searches e + e" — > l + l~~ — > l + l~x 0r X° to derive similar mass limits ( ^80 
GeV) jSH]. A recent quantitative analysis of the above decays in the context of a Little Higgs 
Model can be found in [37} . 

Concerning if ±=t , LEP searched for both left-handed and right-handed (which 
we will not consider in this paper) via several mechanisms: 

(i) Pair production via e + e~ — > j*,Z* — > H ++ H followed by decay to (l^ = 
e , n , r ); the cross-section is determined by gauge couplings and leads to mass limits 
of m H ±± > 100 GeV 0, 0.D33- 

(ii) Single production of if ±± via e + e~ — > H ±± e T e T ; the rate is determined by the cou- 
pling /in and leads to excluded regions in the plane (/in, m#±±), with sensitivity up to 
m H ±± ^ 180 GeV. Limits of 10~ 2 -> 10" 1 were set on h n P]. 

(iii) The effect of if ±± on Bhabha scattering e + e~ — > e + e~; as in (ii) above this leads to 
excluded regions in the plane (hn,m H ±±) [TU], [TT] with sensitivity up to m H ±± <^2 TeV. 
Limits of 10~ 2 — > 10 _1 were set on hn. 

The direct searches for i7 ±± will continue at the hadron colliders, Tevatron and LHC. 

3 Production of H ±i: at the Tevatron 

A distinct signature of if ±=l= would be a pair of same sign charged leptons (e or fi^) with high 
invariant mass. At hadron colliders such a signal has a relatively high detection efficiency and 
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enjoys essentially negligible background from Standard Model processes. Earlier theoretical 
studies of the search potential for if ±± at such colliders can be found in [201,121], with a recent 
analysis at the LHC in j3H] . The decays of H ±± to states involving r are more problematic at 
hadron colliders, although simulations in these channels j2I], [HE] promise sensitivity to values 
of m H ±± beyond the LEP limits. The decays if ±=t — > W ± W ± are proportional to and can 
be neglected in the case of very small v& of interest to us. 

In 2003 the Tevatron performed the first search for if ±± at a hadron collider. DO [IS] have 
searched for if ±=t — > p^pT while CDF searched for 3 final states: if ±± — > e ± e ± , e^/i^, /i ± /x ± . 
The assumed production mechanism for if ±± is qq — > 7*, Z* — > H ++ H . 3 This cross-section 
depends on only one unknown parameter, m#±±, and importantly is not suppressed by any 
small factor such as a Yukawa coupling hij or a triplet VEV. The search assumes that if ±=l= is 
sufficiently long-lived to decay in the detector, which corresponds to hu > 10~ 5 . A search for 
a long lived H ±i: decaying outside the detector has been performed in [39J. The cross-section 
also depends on the hypercharge of the Higgs representation, which is Y = 2 in the HTM. This 
value of Y is also assumed in the experimental searches. The explicit partonic cross-section at 
leading order (LO) is as follows (where q = u,d): 



a LO {qq^H ++ H- 



9Q : 



A 3 



2 2 1 
e q e H + 



e q e H v q v H {l - M 2 Z /Q 2 ) + (v 2 + aj)v| 
(1 - Mf/Q 2 ) 2 + Mfri/Q 4 



(9) 



Here v q = (I 3q - 2e q s 2 v ) / \s w c w ) , & q = h q /{s w c w ), and v H = (I 3H - en-sly) / (s w c w ). The 
third isospin component is denoted by I 3q (hn) and e q (en) is the electric charge of the quark 
q (if ±=t ). sw and Cw are sin 6^ and cos 6w respectively. Q 2 is the partonic centre-of-mass 
(COM) energy, a is the QED coupling evaluated at the scale Q, Mz is the Z boson mass, 
Y z is the Z boson width, and (3\ = \Jl — 4m^ ±± / Q 2 . Order a s QCD corrections modify the 
LO cross-section by a factor K « 1.3 at the Tevatron for m^±± < 200 GeV, and K w 1.25 at 
the LHC for m H ±± < 1000 GeV [10]. We neglect the gluon-gluon fusion (a 2 ) contribution to 
H ++ H production, which has no compensatory enhancement factor analogous to the tan 4 f3 
term for doublet H ± production via gg — ► H + H~ |41j . 

Assuming that if ±± production proceeds via this pair production process, the absence of 
signal enables a limit to be set on the product: 



a(pp -> H ++ H—) x BR(H 



Itl}) 



(10) 



Clearly the strongest constraints on m^±± are obtained assuming BR(if ±=l= — > if if) = 
100%. Currently these mass limits stand at: 133,115,136 GeV for the e ± e ± , e ± p ± 1 // ± // chan- 
nels respectively [TO]. In the HTM one expects BR(H ±i: — > ) 7^ 100% if Eqn.© is required 
to explain the currently favoured form of the neutrino mass matrix |14j . 

The current search strategy is in fact sensitive to any singly produced H ±zL , i.e. signal 
candidates are events with one pair of same sign leptons reconstructing to m#±± . This require- 
ment is sufficient to reduce the SM background to negligible proportions. Hence the search 
potential of the Tevatron merely depends on the signal efficiencies for the signal (currently 

3 The model-dependent contribution from any Z' (which can enhance the cross-section |24 ) .|38 p is currently 
not considered. 
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~ 34%, 34%, 18% for /i/i, ee, e/x) and the integrated luminosity. With these relatively high ef- 
ficiencies and an expected £ = 4 — &fb~ l by the year 2009, discovery with > 5 events will be 
possible for <Jh++h — of a few fb, which corresponds to a mass reach m H ±± < 200 GeV. 

Although single if ±± production processes such as pp — ► W ± — ► W /=F if ±± can be neglected 
4 due to the strong triplet VEV suppression, the mechanism pp — > — > H ±± H T is potentially 
sizeable. This latter process proceeds via a gauge coupling constant and is not suppressed by 
any small factor. The LO partonic cross-section is as follows: 

VLotf* - H++H-) = u ™Iq2 C tPw$ ( n ) 

Here CV arises from the H ±± H T W T vertex and CV = 2 for J = 1,Y = 2 triplet fields (the dou- 
blet component of H ± is negligible) ; (3 2 = ^/ (1 — ijn H ± + m#±±) 2 /C} 2 )(l — (m#± — m H ±±) 2 /Q 2 ) 
and = Q 2 /{Q 2 — Mw)- F° r simplicity, we take the same K = 1.3 as for a{qq — > H ++ H ) 
at the Tevatron and K = 1.25 at the LHC. Explicit calculations for the K factor for the pro- 
cess u{q'q — > H^A ) in the MSSM jH] (which shares the same K factor as q'q — > H ++ H~) 
give K « 1.2. In this paper we will study in detail the magnitude and relative importance of 
a {q'q — > H ±:k H T ). Although we work in the HTM, our numerical analysis is relevant for other 
models which possess a I = 1,Y = 2 Higgs triplet (e.g. L-R symmetric models and Little Higgs 
Models). 

A previous quantitative study of this mechanism can be found in Cross-sections were 
given at both LHC and Tevatron energies for m^±± > 200 GeV with the simplifying assumption 
m H ±± = m H ±. It was shown that a{q'q — > if ±± if T ) can of comparable size to a{qq — > 
H++H-). 

In this paper we first generalize the work of [21] as follows: 

(i) In our discussion at the Tevatron we consider masses in the range 100 GeV < m#±± < 200 
GeV which will be probed during Run II, and allow mild mass splittings |m#±± — m^± | < 
20 GeV. 

(ii) In our discussion at the LHC we consider larger mass splittings |m#±± — m#± | < 80 GeV. 

(iii) For both the Tevatron and LHC we study in detail the relative magnitude of <j{q'q — > 
H^HT) and a{qq -> #++#—). 

Moreover, motivated by the fact that the currently employed Tevatron search strategy is 
sensitive to single production of i? ±± , we advocate the use of the inclusive single production 
cross-section (er#±±) when comparing the experimentally excluded region with the theoretical 
cross-section. This leads to a strengthening of the mass bound for m^±± which now carries a 
dependence on m H ±. We introduce the single production cross-section as follows: 

<Tff±± = cr{pp,pp — > H ++ H ) + a{pp,pp H ++ H~) + a{pp,pp — > H H + ) (12) 
Single production of a right-handed triplet via q'q — > 

W R -> W R H±± S3 and W^W£ fusion g3| can be 

sizeable at the LHC. 
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Figure 1: (a) Single production cross-section of if ±± (cr#±±) at the Tevatron as a function of 
m#±± for different values of m#±. (b) Ratio R as a function of m#±±. We use CTEQ6L1 
parton distribution functions. 

At the Tevatron a(pp -> H ++ B~-) = a(pp -> H~H + ) while at the LHC <j(pp -> H ++ R~-) > 
a(pp — > if If a signal for if ±=t were found in the 2 lepton channel, subsequent searches 

could select signal events with 3 or 4 leptons, in order to disentangle gg — ► H ++ H and 
g'g — > H^H^. In our numerical analysis we utilize the CTEQ6L1 parton distribution functions 
(pdfs) [13]. We take the factorization scale (Q) as the partonic COM energy (y/s). Our results 
for a(qq — > H ++ H ) agree with those in j2l],jinj. Our results for a(q'q — > H ±:k H T ) agree 
with those in [21] (and taking C T = 1 agree with a(q'q -> if ± A°) in the 2HDM/MSSM (Sj). 
The above cross-sections evaluated with MRST02 pdfs agree with those evaluated with 
CTEQ6L1 to within 10% -> 15%. 

In Fig. ^ (a) we plot cth±± as a function of m^±± at the Tevatron for three different values 
of m H ±. We take K = 1.3. The current excluded regions from the e ± e ± , e ± /i ± ,/i ± / u ± searches 
correspond to the area above horizontal lines at roughly 40, 70, 35 fb respectively. The present 
mass limits for m#±± are where the curve for H ++ H intersects with the above horizontal lines, 
and read as 133, 115, 136 GeV respectively for BR( J ff ±± -> if if) = 100% . With the inclusion 
of the i7 ±± i7 =F channel, these mass limits increase to 150, 130, 150 for m#± = m#±± + 20 
GeV, strengthening to 160,140,160 for m#± = m#±± — 20 GeV. Clearly the search potential 
of the Tevatron (i.e. the mass limit on m#±±) increases significantly when one includes the 
contribution to cr H ±± from pp H^H^. Note that the above mass limits strictly apply to the 
case when H ±±L decays leptonically, and with BR=100% in a given channel. However, if hy are 
to provide the currently favoured form of the neutrino mass matrix then BR(if ±=l= — > if if) < 
100% in a given channel. Moreover, if m^±± > m^± then the decay channel H ±zL — > H ± W* 
would be open. As shown in (221 > this decay can be sizeable and thus reduces BR(if ±=t — > if if). 
We will return to these issues in Section 5. 

In Fig. ^ (b) we plot the ratio of cross-sections R at the Tevatron as a function of m#±±, 
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where R is defined as follows: 



_ cr(pp,pp -> ) + a(pp,pp -» if ff+) 

a(pp,pp — >■ H ++ H ) 

The m#±± dependence arises from the phase space functions /?i and /3 2 in EqnsE] and [TT] As 
can be seen, 0.8 < R < 2.2 and thus q'q —> H^H^ contributes significantly to <x#±±. 

In Fig. 2 we plot the analogies of Fig.l for the LHC. In Fig. |21 (a) we plot er#±± for 3 values 
of rriH±± and for larger mass splittings (|m#±± — m#±| < 80 GeV) than in Fig. |21 We take 
K = 1.25. As before, the inclusion of q'q — > H^H^ significantly increases the search potential 
e.g. if sensitivity to <Jh±± = 1 fb is attained, the mass reach extends from mn±± < 600 GeV 
(H ++ H only) to 750 GeV for {m H ± = m H ±± — 80 GeV). Recently [SH] performed a simulation 
of the detection prospects at the LHC for qq — > H ++ H for the cases where 3 and 4 leptons 
are detected. With 100 fb -1 , sensitivity to m H ±± ^800 GeV (3 leptons) and m H ±± <^700 GeV 
(4 leptons) is expected. We are not aware of a simulation for the case where only 2 leptons 
are detected. Presumably even larger values of m H ±± ( ^ 800 GeV) could be probed. In Fig. 
12] (b) we plot R as a function of m H ±±. One can see that R > 1 for the upper two curves for 
all m H ±±, while for the lower curve R > 1 for m#±± > 260 GeV. Note that the dependence of 
R on rriH±± differs from that observed in Fig. Q (b), which can be attributed to the different 
parton luminosity functions at the Tevatron and LHC. 




Figure 2: (a) Single production cross-section of if ±± (<t^±±) at the LHC as a function of m#±± 
for different values of m^±. (b) Ratio R as a function of m^±±. We use CTEQ6L1 parton 
distribution functions. 
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4 Neutrino mass hierarchy and the decay H ±zL —> H ± W* 

The current experimental searches assume that the sole decay mode of if ±± is H ±zSz — > if if 
mediated by the arbitrary Yukawa couplings h^. The decay rate for if ±± — > Z^ 1 is given by: 

r(// ±± ^/^) = 5^|^f (14) 

where 5 = 1(2) for z = j [i ^ j). Clearly r(if ±=t — ► Zf if) depends crucially on the absolute 
value of the h^, although the leptonic BRs are determined by the relative values. In this 
section we consider the impact of the decay mode H ±zL — > H^W* on the BRs of the leptonic 
channels. It has been known for some time that BR(if ±=l= — > H W*) is potentially sizeable 
and a quantitative analysis can be found in [22]. The decay rate for J ff ±± -> tf+W* summing 
over all fermion states for W* — > // excluding the t quark) is given by: 

T( J ff ±± -> iT^*) = 9G|M^m H ±±C|P/(167r 3 ) (15) 

where P is the phase space term (which we calculate by numerical integration) and CV(= 2) 
is from the coupling H^H^W . P depends on the mass difference Am defined by Am = 
m#±± — m#±, and P = for Am = 0. If m#± < m#±± this decay can compete with 
i7 ±± — > Zf ^ since the phase space suppression of the virtual W* is compensated by the gauge 
strength coupling [2T]. Ref.j22] showed that H ±zL — > iP^VK* can dominate over if ±± — > Z?^ 
if Am is sizeable (> 40 GeV) and Ztjj are of order 1CT 3 or less. A large BR(if ±± — > H ± W*) 
would debilitate the if ±=t search potential in the leptonic channel. However, as emphasized in 
[2*3] . observation of if ±=t — > H ± W* together with one or more of the leptonic channels could 
provide information on the absolute values of h^. If only BR(if ±± — > Z, lj ) are measured 
then only the relative values of the can be evaluated. The decay rate for if ±± — ► H^W* is 
theoretically calculable once m,H± and m#±± are known experimentally, and thus it can be used 
as a benchmark decay with which to estimate the total width of if ±=t . It is known that the 
BRs of the leptonic channels depend on which solution to the neutrino mass matrix is realized 
|14j . However, a quantitative analysis of the impact of if ±± — > H^W* in the various allowed 
scenarios is still lacking and will be presented below. We are not aware of any experimental 
simulation of if ±± — > H ± W*. The signature would depend crucially on the decay products of 
H ± , which are are either H ± — > / ± z/ i (driven by hij), or possibly H ± — > H°W*, A°W*. 

We now briefly review relevant results and formulae from neutrino physics. The neutrino 
mass matrix is diagonalized by the MNS (Maki-Nakagawa-Sakata) matrix V MNS [47;. Using 
Eq.fJSJ) one can write the couplings h^ as follows: 

h ij = '^^ V M^s dia 9 (mi, m 2 , m 3 )V^ NS (16) 

Here we take the basis in which the unitary matrix responsible for diagonalizing the charged- 
lepton mass matrix is a unit matrix. The MNS matrix in the standard parametrization is as 
follows: 

C\C 3 S1C3 
r |X> : [ -sica - c lS2 s 3 e lS Cl c 2 - s lS2 s 3 e l5 s 2 c 3 | e^ 2 |. (17) 
sis 2 - cic 2 s 3 e lS -cis 2 - sic 2 s 3 e lS 
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where Sj = sin0j and q = cos0j, 5 is the Dirac phase and (fi and ^2 are the Majorana phases. 

Neutrino oscillation experiments involving solar , atmospheric j3Hl and reactor neutrinos 
jSHl are sensitive to the mass-squared differences and the mixing angles, and give the following 
preferred values: 



Am 2 2 



m 



8.0 x 10" 5 eV 2 



\Am 2 13 \ 



m, — m 



2.1 x lO^eV 2 



sin 2 20i ~ 0.8 , sin 2 26 2 



sin 2 20 3 £0.16. 



(18) 
(19) 



Since the sign of Am 2 3 and the mass of the lightest neutrino are both undetermined at present, 
distinct neutrino mass hierarchy patterns are classified as follows: Normal hierarchy (NH) 
(mi < m2 <C ^3); Inverted hierarchy (IH) (rri2 > m\ ^> 7713), Quasi- degenerate (DG) (mi ~ 
m 2 ~ m 3 ^> 



Am 13 1). From Eq.fjSJ) and Eq. (|T6*j) it can be shown that: 



2 

i J 



(20) 



l :J 



Hence the total leptonic decay width depends on the absolute mass of the neutrinos, and the 
value of Yji m li depends on which solution to the neutrino mass matrix (NH,IH,DG) is real- 
ized. The minimum value of J2i m i is |Am 2 3 | while the maximum is given by the cosmo logical 
constraint. 

In FigJ3]we show contours of BR(iJ ±=t — * H ± W*) in the plane (m H ±±, v a), for three different 
solutions to the neutrino mass matrix. We assume that = for NH (IH) and m\ = 0.2 

eV for DG. We take m H ± = m H ±± — 20 GeV. From EqUHJ all are determined once v a is 
specified. In order to comply with current experimental upper limits on LFV decays of fi^ 
and t^, one can derive the bound v a > 10 eV for NH and IH, and v a > 100 eV for DG. The 
stronger constraint on v& in DG arises because J2i m i i n DG is larger than those in NH and IH. 



(a) 

NH 

















>H i W*) >80% 








60% 








40% 








20% 



(b) 

IH 



BR(H" ->H W ) >«0% 



20% 



+ (GeV) 



(c) 
DG 




L 


am* ->HVl >80Si 








60* 








.to'; 



























Figure 3: Contours of BR(iJ ±± -> H^-W*) in the plane (m H ±±,v A ), for NH (a), IH (b) and 
DG (c). We take m H ± = m H ±± - 20 GeV. 



From FigBJit is clear that BR( J ff ±± iZ^W*) can be sizeable, and approaches 100% for 
larger For a fixed value of v&, one can see that BR(if ±=l= — > H ± W*) is relatively more 
important in NH and IH than in DG. This can be understood from Eq. I2H since DG requires 
heavier neutrinos (and thus larger hij) which in turn reduces BR(if ±± — > H^W*). One can 
consider three distinct scenarios with very different magnitudes for BR(ff ±± — ► H ± W*) and 
BR( J ff ±± -> 1*1*): 
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(i) BR( J ff ±± -> iZ^W*) > BR( J ff ±± -> / ± / ± ): 

In this case the current search strategy (which requires if ±± — > Z Z decay) is ineffective. 
Simulations have not been carried out for the decay i7 ±± — > H^W* although one might 
naively expect sensitivity comparable to that for the decay if ±=t — > r ± r ± , as suggested 

in EH, ESI- 

(ii) BR(i/ ±± -> iZ^W*) « BR(if ±± -> Z ± / ± ): 

The search for if ±=l= — > / ± / ± would be effective and H ±:iz could be discovered in one 
or more leptonic channels. If iZ ±± — > H ± W* is also observed then information on the 
absolute value of ft^- might be possible: Using EqsJ14l and ITol the ratio of leptonic events 
(Nu,) to H ± W* events (N H ± W *) is given as follows: 



N «> (21) 



Observation of the leptonic channel provides m^±±. If m#± can be roughly measured 
then P (and hence the partial width for if ±± — ► H ± W*) can be calculated. From the 
above equation one can obtain an order of magnitude estimate of hij. 

(iii) BR(F ±± -> iZ^W*) < BR(iJ ±± -> / ± Z ± ): 

In this case the current search strategy (H ±i: — > is effective. If BR(iZ ±=t — > Z^/j 1 ) 

are measured then the ratios of Zijj can be evaluated. This can be compared with EqnJ3] 
in order to see which neutrino solution is realized [14] . The absolute values of cannot 
be measured unless a LFV decay of fi and/or r is observed. 



5 Tevatron search potential in HTM 

We now study the search potential of the Tevatron for the generalized case in the HTM where 
pp — > H is included, BR(iZ ±± — > H ± W*) ^ 0% and hij are required to reproduce 

a phenomenologically acceptable neutrino mass matrix. We relax the assumptions for the 
Majorana phases and take tp 2 = or tt, which leads the 7 distinct solutions: 

NH: mi < m 2 <C m 3 , 

IH1: m 2 > mi ^> m^, IH2: —m 2 > mj > m 3 , 

DG1: mi — rri2 — ™>3, DG2: nil ~ rri2 — -m3, 

DG3: mi ~ — m 2 — ms, DG4: mi ~ —m 2 — — m3. 



In the HTM, BR(iZ ±± — > Z ± / ± ) are predicted and different in each of the 7 distinct solutions 
(NH,IH1,IH2,DG1— >DG4), and their ratios were evaluated in Note that such predictions 
of BR(iZ ±± — > / ± / ± ) are a feature of the HTM in which the couplings are the sole origin of 
neutrino mass. This direct correlation between BR(iZ ±:t — > 11) and the neutrino mass matrix 
may not extend to iZ ±=t of other models in which neutrinos can acquire mass by other means 
e.g. the seesaw mechanism in L-R models or by a combination of mechanisms which may or may 
not include the couplings jSZ]j|SIl- In contrast, the production process cr(pp — > H^H^) is 
certainly relevant in any model with Y = 2 triplets. 
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In Figs|U— > IH1 we plot an as a function of m H ±±, where an is the total leptonic (/ = e,//, r) 
cross-section defined as: 

<r w = <r(pp -> H ++ H—) x B«(2 - B«) + 2a(pp -> #++#-) x (22) 

The contribution to cr^ from cr(pp — > H ++ H ) falls more slowly with decreasing Bn since 
signal candidates are events with at least 2 leptons. Eq. (|22|) simplifies to EqEl in the limit 
where a(pp — ► H ±:k H T ) = and S/z = 1. FigsEfa) shows 07; for the NH with m#± = m H ±±, 
which leads to Bu = 1. In this case X)°"zz — Gh±±- F° r the other figures we take m#± = 
m#±± — 20 GeV, which induces a sizeable (but not dominant) BR(if ±± — > H ± W*), and hence 
X)o"zz < cTff±±- We set i>a = 10 eV in Figs0] and 5 and v& = 100 eV in FigsEl We only 
plot an f° r ee i since the Tevatron has already performed searches in these channels. 

Sensitivity to an of a few fb will be possible with the anticipated integrated luminosities of 
4 — 8 fb -1 . There are plans to search for the 3 leptonic decays involving r (ex, fir, rr) although 
the discovery reach in m#±± is expected to be inferior to that for the ee, e/i, /i/i channels. In 
all figures we take 6 3 = 0°. From the figures it is clear that a eefi ^^ differ considerably in each 
of the 7 scenarios. Optimal coverage is for cases DG1 and DG4, which have a ee ^ > 5 fb and 
a en,nii > 5 fb respectively for m H ±± ^ 180 GeV. For NH, a m > 5 fb for m H ±± <; 190 GeV but 
a ee and a efl are both unobservable. Taking # 3 at its largest experimentally allowed value results 
in minor changes to all figures, with the most noticeable effect being a significant reduction of 
a^ in DG4. Clearly the Tevatron Run II not only has strong search potential for H ±zil , but 
is also capable of distinguishing between the various allowed scenarios for the neutrino mass 
matrix. 




Figure 4: an as a function of m#±± for NH with (a) m#± = m#±± and (b) m#± = m#±± — 20 
GeV . 
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Figure 5: an as a function of m H ±± for (a) IHl and (b) IH2. 

6 Conclusions 

We have studied the production of doubly charged Higgs bosons (-ff ±± ) at hadron colliders 
in the Higgs Triplet Model (HTM), in which a complex Y = 2 scalar triplet is added to 
the Standard Model. The HTM can explain the observed neutrino mass matrix by invoking 
Yukawa couplings hij of the triplet fields to leptons. A definitive signal of the HTM would 
be the observation of the decay if ±=t — > / ± / ± , which enjoys almost negligible background at 
hadron colliders, and whose branching ratios are correlated with the neutrino mass matrix. We 
studied the production mechanism q'q — > H ±:iz H T which can be as large as the mechanism 
qq — > H ++ H assumed in the current searches at the Tevatron. Since the present search 
strategy is sensitive to single production of -ff ±=t , we advocated the use of the inclusive single 
production cross-section (cr#±±) when comparing the experimentally excluded region with the 
theoretical cross-section. This leads to a strengthening of the mass bound for m#±± which 
now carries a dependence on m H ±, and significantly improves the if ±=t search potential at the 
Tevatron and LHC Although we performed our numerical analysis in the HTM, we emphasized 
that the introduction of 0"#±± is also relevant for any model which contains a Y = 2 Higgs triplet 
(e.g. L-R symmetric models and Little Higgs Models). 

Moreover, we quantified the impact of the decay mode iif ±± — > H ± W* for the case of a 
hierarchical, inverted hierarchical and degenerate neutrino mass spectrum. On discovering a 
H ±i: it would be imperative to measure the absolute value of (and hence v a) in order to 
reconstruct the low energy Higgs triplet Lagrangian. We stressed that an order of magnitude 
estimate of could be obtained if the channel if ±=t — > H^W* is observed andrriH± is roughly 
measured. We encourage a detailed experimental simulation of this decay mode at both the 
Tevatron and LHC. 
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Figure 6: on as a function of m H ±± for (a) DGl, (b) DG2, (c) DG3 and (d) DG4. 
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